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We view MgBi as electronically equivalent to (non-staggered) graphite (B~ layer) that has 
undergone a zero gap semiconductor to a superconductor phase transition by a large c-axis (chemical) 
pressure due to Mg ++ layers. Further, like the pir bonded planar organic molecules, graphite 
is an old resonating valence bond (RVB) system. The RVB's are the 'preexisting cooper pairs' 
in the 'parental' zero gap semiconducting B~ (graphite) sheets that manifests themselves as a 
superconducting ground state of the transformed metal. Some consequences are pointed out. 
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Nature has once again surprised us through a high 
temperature superconductivity in MgB% , discovered by 
Akimitsu and collaborators JjJ. This has resulted in a 
spurt of activity both in the experimental and the- 
ory [p~p| p~2|| front. The aim of the present letter is to 
search for possible electronic contribution to the mech- 
anism of superconductivity in MgB 2 ■ Our theory uses 
the well recognized connection of electronic structure of 
boron sheets to graphite sheets ppc|]. In particular we 
view MgB~2 as a (non-staggered) graphite that has un- 
dergone a semi-metal to superconductor transformation 
because of a c-axis chemical pressure. 

The following are some of the phcnomcnological rea- 
sons for looking for an electronic contribution to super- 
conductivity in MgBi : i) similarity to pn bonded pla- 
nar systems like graphite, organic molecules and carbon 
nano-tubes where correlations play varying roles ii) ab- 
sence of Hebel- Slichter peak in NMR relaxation || iii) 
a temperature dependent peak around 17 meV in en- 
ergy resolved neutron scattering B iv) first order metal 
to metal transition on Al or C substitution j7| v) pos- 
sible anomalous temperature dependence of Rh |)j and 
also London penetration depth || and vi) need to use an 
anomalously small fi* in phonon based theories |r| to 
get a reasonable T c . 

MgB~2 is very similar to graphite both crystallographi- 
cally and electronically. While the carbon hexagonal lat- 
tices are staggered in graphite, the boron hexagonal lay- 
ers are on top of each other in MgB 2 . The boron layers 
alternate with a triangular lattice of Mg layers. There is 
nearly complete 2e charge transfer from Mg to the boron 
subsystem: MgB 2 = Mg ++ (B~) 2 - Each boron acquires 
one electron and acquires the electron configuration of a 
carbon atom: 

B-(2s 2 2p 2 ) = C(2s 2 2p 2 ) 

Thus the B~ sheets are electronically like graphite sheets. 
The Mg ++ ion with its strong positive charge pulls the 
charged boron sheets closer and reduces the c-axis sepa- 
ration. The ratio of B-B distance along the ab plane to 
that along the c-axis is w 2.0 as opposed to the corre- 
sponding C-C distance ratio « 2.4 in graphite. This 15 



% shortening of c-axis distance and removal of stagger- 
ing converts a serai conducting 'graphite' into a high Tc 
superconductor ! 

Graphite, like the pir bonded planar organic molecules, 
has been the testing ground for the ideas of resonating 
valence bonds (RVB) from early times. While Pauling 
fuj and others emphasized the idea of valence bond res- 
onance and derealization in the calculation of ground 
state (cohesive) energy and to some extent a quantitative 
understanding of anomalous diamagnetism |l5| ] , non triv- 
ial consequences of RVB ideas for graphite have not been 
seriously discussed. The aim of the present letter is to ar- 
gue that RVB's, the preexisting enhanced singlet (cooper 
pair) correlations in graphite, a zero gap semiconductor 
should reveal itself as a superconducting ground state af- 
ter an application of sufficient c-axis pressure (when no 
structural modification intervenes). Nature seems to re- 
alize this through a chemical pressure in MgB 2 , a system 
isoelectronic to graphite. 

In order to understand MgB 2 we should understand 
an isolated graphite layer. The relevant valence orbital 
is the p z orbitals of carbon with a mean occupancy of 
one electron. In the honeycomb structure there are two 
equivalent carbon atoms and two valence electrons per 
unit cell. According to band theory, in view of two or- 
bitals per unit cell, we get an empty and a filled band. 
Symmetry forces either a zero or a positive overlap be- 
tween the two bands. That is we either have a zero gap 
semiconductor or a semi metal. For a pir bonded sys- 
tem like graphite a well known model is the Hubbard 
model (simplified PPP model) on a honeycomb lattice 
with nearest neighbor hopping: 



H 



-t^ClC^ + h.c. + U 

(ij) 
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The kinetic energy term gives us a zero gap semiconduc- 
tor with a valence and conduction band dispersion 



e±(fc) =±2t 
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of the C layers with t w 2.5eV. The zero gap is not 
due to electron correlation; it is symmetry dictated. It 
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disappears and produces small electron and hole fermi 
pockets if we introduce next nearest neighbor hopping 
term in the planes or finite c-axis hopping term t± . Ex- 
perimentally, in pure graphite the electron and hole con- 
centrations are rather small ~ 10 _4 electron or hole per 
site; pure graphite sheets cleave and the binding between 
layers is van der Waals and not metallic. 

It is known that in pn bonded planar organic 
molecules the screened on site coulomb interaction in car- 
bon (and also boron) p z orbitals is U ~ 6 eV (the bare 
atomic U ~ 12 eV), making the ratio y ~ 2.5. Accord- 
ing to early studies of Sorella and Tosatti and a very 
recent study of Furukawa 0] a Mott insulating behavior 
is obtained for a honeycomb lattice only when — > 3.5. 
Below this critical value of -j the zero gap character of 
the one electron states should make the many body ef- 
fects at low energies less pronounced even though one is 
in 2 dimensions. However, some studies fl7| , |l8| , indicate 
anomalous life time for the quasi particles close to the 
fermi level in graphite. 

The intermediate coupling character of the Hubbard 
model for graphite makes it very difficult to approach it 
either as a weak coupling or a strong coupling problem 
analytically. To circumvent this problem, and also in- 
spired by early RVB ideas on graphite, we propose an 
effective Hamiltonian for graphite/Mg^ on semi phe- 
nomenological grounds, and later sketch a possible mi- 
croscopic derivation of this from the Hubbard model. 

In the early treatment of graphite and pir bonded pla- 
nar organic molecules, such as the discussion due to 
Pauling |l4|], the resonating valence bond (RVB) and en- 
hanced valence bond amplitude is emphasized. Configu- 
rations with nearest neighbor singlet bonds (VB) are en- 
couraged in comparison to the polar (double and single 
occupancy in the p z orbitals) configurations. This con- 
cept gave good estimate for cohesive energy, C-C bond 
distance and even some excited state properties, such 
as the singlet triplet exciton energy differences. From 
the theoretical study |]l9|-^l| of cuprates in the last one 
dozen years we have also realized that RVB's are also 
preformed cooper pairs. They may have profound conse- 
quences such as high temperature superconductivity un- 
der appropriate conditions. 

Our primary aim now is to incorporate the well recog- 
nized RVB character of graphite in our effective Hamil- 
tonian for MgB2 ■ With this in mind we propose the 
following model Hamiltonian for MgBi : 

i,n,a 

~J ^in,jnPin,jn ^ ^-L ^""1 ^in.in+l^in,in+l (^) 

(ij),n i,n 

Here i or j denotes a lattice point on the 2d honeycomb 
lattice and 'n' is the layer index. The singlet operator 
fyj = "^(CjfCji — The two body term with 



Jij > represents an energy gain when two electrons in 
neighboring sites i and j form a singlet (valence bond), as 
is a number operator for bond singlets. That is, this 
term stabilizes covalent configurations relative to ionic 
configurations between any two neighboring sites. Tight 
binding fit of the states close to the fermi level as obtained 
from ab-initio band structure calculation for MgB2 give 
the values of the parameter t w 1.6eV, t± w 1.25eV. 

Now we outline an approximate microscopic derivation 
of the interaction term of our Model Hamiltonian (equa- 
tion 3) and also estimate the phenomenological param- 
eters J and Jj_. In spirit it is similar to superexchangc 
theory, but it is also designed to handle the weak and 
intermediate coupling regions approximately. 

In graphite as well as MgB^ we are exactly at half fill- 
ing in the sense of having one electron per p z orbital. It is 
here we expect maximum effects of two body interactions 
(Hubbard U) in the ground state. In a free fermi gas in 
tight binding system, the spin state of two electrons in 
neighboring sites can be a singlet or a triplet. The local 
singlet character in k-space, as enforced by Pauli princi- 
ple, does not enforce any type of singlet correlations in 
real space. However, two body collisions arising from the 
Hubbard U term, encourages singlet amplitude between 
two electrons of neighboring sites. This is the origin of 
the kinetic or superexchange term in a Mott insulator at 
and close to half filling. 

As in superexchange perturbation theory, we concen- 
trate on two sites involving two electrons, but solve the 
two site problem exactly. The two electron ground state 

is a singlet with an energy E g = — h [U 2 + 16i 2 ] 2 + ¥. 
The energy of the triplet state is at Et = 0. Two ex- 
cited singlet states are at positive higher energies. In 
our effective Hamiltonian (equation 3) the energy gain of 
a nearest neighbor singlet (valence or covalent) bond is 
represented by the term 



Jbtb. 



where J = - [U 2 + 16t 2 ] ' 



U 
2~ 



(4) 



In the large t >> U limit the singlet stabilization energy 
J s» -jj-, agreeing with the superexchange perturbation 
theory. 

For MgB2 using equation (4) we estimate J ~ 1.3eV 
and J_l ~ OAeV. Ideally, in the conducting metallic 
state it is more meaningful to talk about residual cou- 
pling among excitations in momentum space. From this 
point of view the above estimates of singlet stabilization 
terms in real space should be treated with caution. 

The valence bond stabilizing term is an attraction in 
the two body BCS singlet channel, near the fermi surface 
for our half filled case: 



Jj>t 



' J 2li C k^-kl C -k'i C k"[ 



(5) 
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This attraction, under normal conditions, will lead to a 
superconducting ground state. However, graphite is un- 
usual in the sense the density of states vanishes at the 
fermi level. A simple cooper pair analysis shows us that 
the above attraction is incapable of converting the zero 
gap semiconducting state into a superconducting state 
unless the attraction J exceeds a critical value J c w 3t. 
Thus graphite does not have a superconducting ground 
state, in spite of its RVB character ! 
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FIG. 1. Temperature - (c-axis) chemical pressure phase 

diagram. Dashed line is the first order phase transition line 

that ends in a critical point 

Now we qualitatively discuss the zero gap semicon- 
ductor to superconductor transition when we go from 
graphite to MgB 2 ■ As mentioned earlier, the ratio of 
the B-B distance along the c-axis to that in the plane 
gets reduced from 2.4 to 2.0, when we go from graphite 
to MgB 2 . The coupling along the c-axis is van der Waals 
like in graphite. Early band structure calculations and 
transport studies have suggested a c-axis hopping matrix 
element t± < 0.25eV for graphite. The band structure 
results for MgB 2 show a large value of t± <~ 1.2eV re- 
sulting in two large pn fermi surfaces and an addition of 
a small hole density to the the 2 dimensional a bands of 
the honeycomb lattice. 

Thus the major difference between graphite and 
MgB 2 is the c-axis metallization resulting from the chem- 
ical pressure. This has given rise to a finite density of 
states at the chemical potential for MgB 2 ■ Further, the 
local electronic structure and lattice structure has not 
undergone any qualitative modification between graphite 
and MgB 2 ■ Thus MgB 2 , in view of its finite density of 
states at the fermi level is capable of utilizing the local 
RVB correlations and go into a superconducting state. 

The vanishing density of states due to the zero gap 
character in MgB 2 makes the valence bond stabilization 
term 'irrelevant', as long as J < J c « 3t. That is, as a 
function of J there is a quantum phase transition from a 
zero gap semiconductor to a superconductor. The value 
we have estimated for graphite J << 3i, is consistent 



with the experimental observation that graphite is not a 
superconductor. 

A simple BCS mean field theory of our model Hamil- 
tonian (equation 3) gives the approximate formula for 
T c 

k B T c KWe~^k , (6) 

Where W is the bandwidth of the pir band and p is 
the density of state at the fermi level. From the band 
structure calculations, there are two pir bands of width 
18 and 13 eV. We can take the mean bandwidth to be 
W « 15eV~. The mean density of states arising from the 
two bands is po ~ ^ states per eV. If we assume a value 
of J w 0.9eV we get a T c AQK. This suggests that our 
renormalized J s=s 1.3, as obtained from our microscopic 
derivation is in the right ball park. 

In our model, as we increase t± continuously from zero, 
the T c increases continuously; in reality, with the increase 
of c-axis pressure we expect a first order phase transi- 
tion. The zero gap semiconducting state has a screening 
which is fundamentally different from the metallic screen- 
ing of a semi metal with a finite band overlap. On general 
grounds, using arguments similar to Mott, the metaliza- 
tion along the c-axis arising from an uniaxial pressure will 
be a first order phase transition. The c-axis lattice pa- 
rameter collapse arising from the c-axis metalization will 
add to this and make it a stronger first order transition 
like some of the known metal insulator transition. 

We have sketched our general view on the zero gap 
semiconductor to superconductor transition schemati- 
cally in figure 1 and 2. In figure 1, there is a first order 
zero gap (or very small gap) semiconductor to supercon- 
ductor phase transition as a function of pressure in the 
ground state. 
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FIG. 2. Chemical potential (/i) - (c-axis) chemical 
pressure phase diagram. The line separating graphite and 
MgB2 is a first order phase boundary and it ends at a critical 
point 

With a small addition of carriers (through Al substi- 
tution at Mg site for example) we expect the first order 
phase transition to survive as a function of chemical pres- 
sure. In figure 2, we have extended the phase diagram 
to the chemical potential (related to doping) vs chemical 
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pressure (related to the c-axis lattice parameter) phase 
diagram. The first order line ends in a critical point. 

Experiments JjJ on Mg\_ x Al x Bi and MgB2- x C x have 
been performed for a range of x. In the case of Al sub- 
stitution, T c decreases gradually until x 0.1 and for 
0.1 < x < 0.25 we have a two phase region. Supercon- 
ductivity disappears across the two phase region and the 
c-axis lattice parameter jumps down by about 10 %. Al 
being trivalent adds an additional electron to the Boron 
subsystem. We believe that with Al or C substitution 
one crosses the first order phase boundary. 

It is also known that intercalated graphites show low 
temperature superconductivity. According to our picture 
intercalated graphites never get a large fermi surface (and 
hence low T c ) ; as the average c-axis separation increases 
on intercalation, the c-axis metalization does not occur. 

There are several issues that need to be understood 
both qualitatively and quantitatively from our stand 
point. Further our theory has certain consequences which 
can be tested experimentally. 

Symmetry of the order parameter is an important is- 
sue. Even in cuprates, from the numerical analysis it 
is clear that the d-wave and extended-s wave state are 
nearly degenerate. In the present situation, where dou- 
ble occupancy is not be completely projected in the low 
energy subspace, both s-wave and d-wave solutions are 
possible. The absence of Hebel-Slichter peak in the NMR 
does indicate some similarity to the cuprates and or- 
ganic superconductors. The penetration depth study 
Panagopoulos and collaborators || also points to some 
gapless region. Further studies are necessary to settle 
this issue. 

There has been an intriguing neutron scattering result 
[^j, where one sees two unexpected peaks at 17 meV and 
30 meV in the phonon density of states. The 17 meV 
peak in particular has a temperature dependence that 
peaks around the T c . As the authors have mentioned, 
is there a connection of this with the 41 meV peak of 
YBCO materials ? We suggest that this peak is mag- 
netic in origin and expresses an underlying RVB char- 
acter, and polarized neutron scattering experiment can 
prove its spin 1 character. It is also interesting that this 
peak starts building up around 150 K, indicating some 
kind of preformed pairing activity. It is also possible that 
the intrinsic T c of MgB 2 is larger « 200K and that the 
small carrier concentration in the 2d a bond is acting 
like a source of dissipation for superconductivity in the 
pir system, through interband scattering thereby reduc- 
ing T c considerably. For us the lightly doped 2d a band, 
that is predicted by band theory and used by electron 
and phonon based theories, does not play a crucial role 
in establishing high temperature superconductivity. As 



mentioned earlier, it might interferes with superconduc- 
tivity in the pir band. 

It is also interesting that the Hall resistivity Rh has 
a temperature dependence || below about 150 K. The 
strong temperature dependence of Rh in cuprates are 
known to arise from the built up of RVB singlet cor- 
relations, suggesting a built up of RVB correlations in 
MgB 2 . 

I thank R. Shankar (Madras) for discussions. 
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